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Abstract 

The Liouville equation of a two-level atom coupled to a degenerate bimodal 

lossy cavity is unitarily and exactly reduced to two uncoupled Liouville equa- 

tions. The first one describes a dissipative Jaynes-Cummings model and the 

other one a damped harmònic oscillator. Advantages related to the reduction 

method are discussed. 
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I. INTRODUCTION 



Over the last thirty years many theoretical studies have been done in order to under- 
stand the dynamical nonlinear behavior of an atom in a high-Q cavity. The interest toward 
this research area has been mainly spurred by the huge number of experiments revealing 
interesting features of the quantum radiation-matter coupling [1]. Many theoretical and 
experimental activities aimed, in particular, at understanding simple non trivial models of 
quantum òptics involving a single few-level atom and one or more near resonant modes of the 
quantized electromagnètic field of a cavity. The prototype of such models was that proposed 
by Jaynes and Cummings (JCM) in 1963 [2] and describes a two-level atom resonantly in- 
teracting with a single-mode field. The success of this exactly solvable JCM has stimulated 
an intense research devoted at highlighting and generalizing the original idea and physical 
scenario [3]. 

In this paper we concentrate our attention on the coupling between two degenerate 
mode and an effective two-level atom. The dynamics of this extended JCM model, under 
the hypothesis of ideal cavity, has been recently studied [4] bringing into light the existence 
of purely and attractive quantum effects. It should be noticed, however, that in realistic 
situations, cavities currently used possess a finite quality factor Q so that it is of interest 
to know whether and how the predictions of the aforementioned theory are affected by the 
relaxation of photons in the cavity. In order to do this it is very usual to suppose the system 
coupled to a reservoir represented as a bath of harmònic oscillators, describing its dynamics 
by means of a master equation for its density operator. 

In this paper we will show that the resolution of this problem may be exactly traced back 



3 



to that of two simpler decoupled physical fictitious subsystems. We will indeed demonstrate 
that the master equation describing the dynamics of a two-level atom interacting with a 
degenerate bimodal lossy cavity is unitarily equivalent to a system of two uncoupled and 
solvable Liouville equations. The first one governs the dynamical behavior of a dissipative 
one-mode JC model whereas the other one characterizes the dynamics of a damped harmònic 
oscillator. As we shall see, this kind of reduction of the original mathematical problem 
provides a systematic and simple way to express the mean value of a genèric operator of the 
tripartite system under scrutiny, in terms of expectation vàlues relative to physical quantities 
formally pertaining to the two fictitious subsystems. 

II. THE MODEL AND ITS DYNAMICS 

The model we consider consists of two electromagnètic modes of a cavity interacting 
with an effective two-level atom of transition frequency üu . The two independent cavity 
modes possess the same frequency cu ~ uç, but they differ by polarization or direction of 
propagation. The effective two-mode JC Hamiltonian model describing, in the Rotating 
Wave Approximation (RWA), such a tripartite system can be cast in the following form 
(ft = l): 

2 2 

H = J2 w «í a M + ^oSz + X] 9n {al s - + a i* s +) » (!) 

fl=l i=l 

In equation (1) the atòmic degrees of freedom are, as usual, represented by Pauli pseudo spin 
\ operators while a M (ajj is the annihilation (creation) operator of the /i-th cavity mode. 
The two real constants g\ and g 2 measure the intensity of coupling between the atom and 
the cavity mode denoted by 1 and 2 respectively. 
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Let's now observe that the cavities at disposal of experimentalists are characterized 
by a finite quality factor Q. Realistic vàlues of Q have order of magnitude going from 
Q = 10 8 , in correspondence to operi Fabry-Perot resonator [5] to Q = 3 x 10 10 for closed 
cavities [6]. In order to take into account from the very beginning the presence of cavity 
damping, we may suppose that our system is coupled to a reservoir described as a bath of 
harmònic oscillators [7]. We moreover assume, as usual, that the hypotheses under which 
the Born-Markov approximation may be adopted, are satisfied [7]. Under these conditions, 
the master equation for the density matrix p of the combined atom - two mode field system 
can be written in the form [8,9] 

^ = -i [H, p] + Cp, (2) 

where H is the Hamiltonian given by eq. (1) whereas, assuming zero temperature, the losses 
of the cavity are represented by the superoperator £ defined as 

2 

£p = k Y. ( 2c VP°Í - a l a »P - P a l a i) (3) 
n=i 

where k is the damping constant. In eq. (3) we have realistically imposed that the damping 
constants relative to the two cavity modes are equal. 

III. UNITARY DECOUPLING OF THE LIOUVILLE EQUATION 

In order to solve equation (2), we exploit the unitary transformation realized by the 
operator U defined as 

JJ _ e ~l(al ai -a\a 2 ) ^ 

with 7 free real parameter to be appropriately fixed. 



It is possible to demonstrate [4] that the operator U has the following properties: 

U^aïU = cos + siri7a 2 

U^a 2 U = — siii7ai + cos7a2 (5) 

Taking into account equation (5), it is easy to prové that, transforming the Hamiltonian 
H by means the unitary operator U defined by equation (4), one obtains 

H= U ] HU = H 1 + H 2 , (6) 

Hi= coal^ + lü S z + g eff (a\S- + aiSV) ( 7 ) 

H 2 = uala 2 , (8) 

where g eff = \J gj + g\, provided that 7 = tan" 1 (g 2 /gi)- 

The transformed Hamiltonian H describes two independent subsystems in the sense 



that 



Hi,H 2 = 0. Hi describes a simple (single mode)-(single atom) system consisting 



of a fictitious radiation mode, represented by the boson operators a\, a±, linearly coupled 



to the old two-level atom with an effective coupling constant g e ff = y gj + g 2 - H 2 , on 
the other hand, describes a new fictitious radiation mode, with creation and annihilation 
operators given by a 2 and a 2 respectively, decoupled both from the atom and from the first 
collective mode. In other words the canonical transformation of H accomplished by U clearly 
brings into light the fact that the two-level atom induces mode-mode coherence properties 
responsible for the collective behavior of the field subsystem of eq. (1). 

The circumstance that the Hamiltonian model (1) is unitarily equivalent to an one- mode 
JCM plus a free collective mode, turns out to be the key for solving the master equation 
(2). In order to better clarify this point let's start by observing that, if A denotes the 

6 



operator obtained transforming by U a genèric operator A of the system under scrutiny, 
that is A = U^AU, the master equation (2) can be equivalently written in the form: 

+ WCpU (9) 

Exploiting eq. (3) and (5), it is possible to demonstrate that 

2 

U*CpU = fc £ (2U f a^UpU^alU - XJ^a\aJÜ~p - pU^a^U) = dp + C 2 p (10) 

where 

= k {2aipa\ - a\aip - pa\a^ (i = 1, 2) (11) 

Suppose now that the initial conditions imposed to the system are such that the trans- 
formed density operator p(0) can be factorized into a product of two contributions, p~l relative 
to the atom and the collective mode 1 and p\ describing the free fictitious mode 2. In ref. 
[4], for example, the authors suppose that at t — the system is prepared putting the atom 
in its ground state, exciting one mode in a coherent state \a) and leaving the other one in 
its vacuum state. It is easy to demonstrate that such an initial condition, in correspondence 
to which the dynamics of the system manifests new and interesting nonclassical features, is 
such that the condition p(0) = pjpg is satisfied. Generally speaking, however, even if the 
initial condition imposed to the system is such that the three subsystems, atom, mode 1 and 
mode 2, are factorized, the transformed density operator p(0) at t — could not maintain 
the same factorization property. 

In the context of this paper we assume that 

p(o) = pIpI (12) 
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H, p 



Our expectation, both on mathematical and physical grounds, is that the solution of eq. 
(9) satisfying such a factorized initial condition exists and is unique. 

In a lossless situation eqs. (6)- (8) guarantee that the density matrix operator at a genèric 
time instant t keeps its initial factorized form. We may wonder whether the solution of 
equation (9) satisfying the initial condition given by eq. (12) still possesses such a property. 
We may indeed show that this is the case. To this end we state and prové the following 
Theorem: The unique solution of equation (9) satisfying the initial condition (12) may be 
given in the form p(t) = Pi(í)p 2 (í) where Pi(í) and p 2 (í) are solutions of the two following 
Cauchy problems: 

dpi 



dt 

dp 2 
dt 



—i 



H u pij +dpi, pi(0) = pi 

H 2 , p2 + £-2p2 P 2 (0) = p 2 



(13) 
(14) 



Proof. Looking for solutions of equation (9) in the factorized form p(t) = pi(í)p 2 (í) yields 
[-^ + i[H 1 , pi] - £ipi]p 2 + Pi[-0jr + * [^2, p 2 ] - £ 2 p 2 ] = (15) 



As first step we now prové that eq. (15) necessarily implies 

dpi 



dt 
dp 2 
dt 



+ i 
+ i 



Hi, pij - Api = Aipi, 

H2, p2 — ^2p2 = A 2 p 2 



(16) 
(17) 



with Ai = — À 2 



To this aim, for the sake of simplicity, we rewrite equation (15) in the following form: 



A X B 2 + A 2 B 1 = 



(18) 



where the operators Ai = + i Hi, pi — CíPí and Bi = pj (i — 1, 2) are defined in the 



Hilbert space of the system i. We now prové that, for each fixed i {i = 1,2), it must be: 



Ai = \ t B t (19) 

To this end let's consider a genèric basis {\rii) ® \n 2 )} of the Hilbert space of the total system 
obtained as a tensorial product of two basis, and {(^2)} relative to the Hilbert space 

of the system 1 and 2 respectively. The action of both members of equation (18) on a genèric 
state vector \rii) <g> \n 2 ) of the prefixed basis in the Hilbert space of the total system (l) + (2) 
may be formally put as follows: 

(A 1 B 2 )(\n 1 ) <g> \n 2 )) + (A^dm) <g> \n 2 )) = 

= (AilnOXflalna» + (A 2 |n 2 ))(fli|m)) = (20) 

Suppose now A\ 7^ \\B\, in contradiction with eq. (19). Then there exists at least one 
vector, say \fi\) such that 

A^n,) ^ X.ln,) (21) 

Accordingly equation (20) may be cast in the form 

(A^n^B^)) + (A 2 \n 2 ))(B 1 \n 1 )) = (22) 

whatever \n 2 ) is. This implies that for any \p 2 ) 

(p 2 \B 2 \n 2 ) = (p 2 \A 2 \n 2 ) = (23) 

which in turn leads to the conclusion A 2 = B 2 = evidently false since Tr{B 2 } = 1 at any t. 
Thus Ai = À1-B1 and, repeating the same argument after exchanging 1 with 2, A 2 = \ 2 B 2 . 
It is moreover easy to verify with the help of eq. (20), that necessarily À x = —\ 2 = À. This 
concludes the first step of our proof. 
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The second and last step in the proof of our initial statement consists in proving that 
the only physically admissible solution is that correspondent to À = 0. Let's indeed consider 
the general equation 

^ + *[H, p]-Cp = \p (24) 

Taking the trace of both members of this equation we obtain: 

^Tr{p} - Tr{Cp} = \Tr{p} (25) 

since Tr{[H, p]} = 0. Exploiting, moreover, elementary algebraic properties of the trace 
and taking into account the explicit form of the superoperator C, given by eq. (3) ignoring 
the subscript i, it is possible to demonstrate that also Tr{Cp} is equal to zero. In view of 
these considerations eq. (25) becomes: 

^Tr{p} = \Tr{p} (26) 

We thus may state that, being Tr{p} = 1, eq. (26) necessarily implies À = 0. 

Before concluding this section we wish to underline that the unitary reduction of the 
master equation (2) into a system of two uncoupled Liouville equations (see eqs. (13) and 
(14)), pavès the way to the exact treatment of the dynamics of the physical system under 
scrutiny. 

Let's indeed consider a genèric operator O = f(ai, Cl|, Oj2i C^2i 

S) of the system. In order 

to calculate the mean value (O) of this operator at a genèric time instant t, we observe that 

(O) = Tr{pO} = Tr{{U^pU){U ] OU)} (27) 

But 
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rfpU = p = p ± p 2 (28) 

and 

U ] OU = U ] f{a u a\, a 2 , 4, S)U = f(U^ ai U, U ] a\U, U ] a 2 U, U ] a ] 2 U, S) = f(a u a{, a 2 , 4, S) 

(29) 

where / is a new function of the operators ai, a\, a 2 , a\ and S now related to the two 
fictitious subsystems 1 (atom - mode 1) and 2 (mode 2) obtained after the transformation 
defined by the unitary operator U . 

Let's now observe that in many physical problems of interest, / has the form of a linear 
superposition of products of operators a iy a\, (i = 1, 2) and S that is 

f(a h a], S)=EE fk9Í(oi, al)gl(S) (30) 

fc i 

where g l {i = 1,2) and g 3 are prefixed functions. Replacing U^OU appearing in equation 
(27) by / as expressed by eq. (30) immediately yields 

(O) ^fkTrip^a^^yTri^^ai)} (31) 
k 

Eq. (31) says that the mean value of the operator O can be expressed as combination of 
mean vàlues of operators relative to the two fictitious subsystems 1 and 2 respectively. In 
any case, whatever the form of the function / is, exploiting the already known solutions p\ 
and p~2, describing the dissipative JC model and the damped harmònic oscillator respectively, 
eq. (31) provides a simpler way to evaluate the mean value of a genèric operator O of our 
original physical system. 
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IV. CONCLUSIVE REMARKS 



In this paper we have considered an effective two-level atom resonantly coupled to a 
degenerate bimodal lossy cavity. Our main and novel result is that, under appropriate 
but general enough initial conditions, the Liouville equation of our tripartite system can be 
unitarily and exactly reduced to a pair of uncoupled and simpler master equations describing 
a damped harmònic oscillator and a dissipative one-mode JC model. 

The success of our procedure is intimately related to the knowledge of a unitary operator 
accomplishing the decoupling of one bosonic degree of freedom from the other dynamical 
variables. 

We emphasize that since both these new fictitious problems may be analytically solved, 
our recipe provides the key for an exact analytical treatment of the original Liouville equa- 
tion. A complete investigation of this dynamical problem is behind the scope of this paper 
and will be presented elsewhere. 

An effective way to simplify the evaluation of expectation vàlues of physical quantities 
of interest, with the help of our method, has also been presented and briefly discussed. 

We conclude underlining that the reduction method reported in this paper might also 
be useful to treat master equations relative to other physical systems. As an example, it 
is easy to verify that our procedure is successfully applicable to the hamiltonian model (1) 
revised including the counter rotating terms. 
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